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ABSTRACT 
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One  quick  and  easy  way  of  estimating  scattering  coefficient  is  to 
assume  that  the  antenna  beam  is  confined  in  a  finite  cone  with  constant 
gain  inside  and  the  range  to  the  antenna  is  constant  for  all  the  ceils 
inside  the  illuminated  area.  Also  the  variation  of  the  0°  is  assumed  to 
be  small  enough  inside  the  illuminated  area  so  one  can  estimate  the  o° 
with  only  the  power  return  measurements. 

The  above  may  not  be  true  depending  on  the  antenna  gain  shapes  or  the 
target  characteristic.  In  this  report,  the  so-called  narrow-beam  approxi¬ 
mation  is  examined  using  a  little  more  realistic  antenna  gain  functions  and 
range  variation  and  also  a  specifica0  variation  function.  The  numerical 
integration  shows  that  the  narrow  beam  approximation  can  give  as  much  as 

r 

3.6  dB  error  for  normal  incidence  to  the  target  with  rapidly  varying^0 
when  the  beamwidths  are  9^^ide.  For  the  incidence  angles  of  10^°<ir  larger 
the  error  is  in  the  order  of  a  dB  or  less  depending  on  the  beamwidths  and 
target  characteristics. 
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1.0  BACKGROUND 


Accurate  estimation  of  scattering  coefficient  for  an  area-extensive  tar¬ 
get  using  the  average  return  power  is  not  a  simple  task  because  it  involves 
an  integral  equation.  Hence,  it  is  usually  simplified  by  using  several 
approximations  and  assumptions.  The  purpose  of  this  report  is  to  investigate 
the  validity  of  these  commonly  used  approximations  and  to  evaluate  the 
possible  error  when  these  simplifying  assumptions  are  used  to  calculate  the 
scattering  coefficient,  for  example  antenna  beamwidths  and  target  scattering 
coefficient  angular  variations. 

Axline  [1]  developed  a  matrix  inversion  technique  to  invert  the  integral 
equation  for  various  cases  of  antenna  gain  shapes  and  o°  variation  functions. 
This  technique  requires  power  return  measured  as  a  function  of  incidence 
angle  at  a  small  angular  interval  and  is  not  readily  applicable  when  this 
is  not  true. 

Stiles  et  al.  [2]  described  an  empirical  procedure  to  estimate  the 
error,  when  the  narrow-beam  approximation  is  used  for  small  incidence  angles, 
by  initially  hypothesizing  the  theoretical  o°(0)  curve. 

Moore  [3]  suggested  an  iteration  scheme  to  estimate  the  scattering 
coefficient  by  using  the  constant-a°  approximation  to  get  a  first  estimate 
of  the  o°  variation  function  and  then  continuously  refining  that  function 
to  get  better  estimate  of  a0.  This  method  is  used  in  this  report,  and  the 
integral  is  evaluated  numerically  using  the  real  system  parameters  of  the 
HELOSCAT  III  FM-CW  scatterometer  and  a  data  set  taken  from  sea  ice. 


2.0  NARROW-BEAM  APPROXIMATION 


Consider  the  radar  equation 


J  J 

area 

i 1 iuminated 


P  g2xV 

— TT  dA 

(W)V 


(l) 


where  identical  transmit-  and  recei ve-antenna  gains  are  used. 

The  first  problem  here  is  how  to  define  the  illuminated  area.  Actually 
the  antennae  have  lobes  everywhere,  so  the  area  has  to  be  from  -®  to  +°°  in 
the  plane  of  illumination.  Usual  practice  for  this  is  to  assume  that  the  an¬ 
tenna  beam  is  confined  in  a  finite  cone  within  which  the  maximum  value  of  gain 
applies,  with  zero  gain  elsewhere.  The  effective  beamwidth,  B  is  defined 
in  such  a  way  that  the  pattern  solid  angles  for  the  two-way  real  antenna 
beam  and  the  cone-shaped  beam  are  the  same  (see  Figure  1).  Therefore, 
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where  Gg  is  the  maximum  gain  for  the  two-way  pattern.  This  is  not  a  simple 
2 

task  unless  G  is  completely  known.  If  the  gain  function  is  assumed  to 


be  a  Gaussian  and  if  it  is  $- independent,  then 


which  reduces  to 


*eff  =  2  cos"1(I 


■2.7736/6 
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This  equation  can  be  solved  numerically  and  for  the  of  up  to  15°,  the 
following  relation  can  be  used. 


eff 


1.2  6 


3dB 
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Antenna  beamwidths  are  usually  measured  in  azimuth  and  elevation  planes. 
When  the  measured  patterns  are  considerably  different  from  each  other,  separate 
effective  beamwidths  are  defined. 


Ba,eff  "  1,2  6a,3dB 
Be,eff  "  1,2  6e,3dB 

The  illuminated  area  is  then  the  ground  projection 
beam  and  the  area  of  the  ellipse  is  approximately 


(6) 

(7) 

(Figure  2)  of  the  conical 
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If  the  beamwidths  are  narrow  enough,  RH  inside  the  integral  in  Eq.  (1) 

can  be  treated  to  be  constant;  and  if  a°  varies  slowly  inside  the  illuminated 

area  given  by  Eq.  (8),  one  can  remove  the  o°  from  the  integral  and  solve  for 
o 

a  . 


o  _  (Att)3  /Pr>  ,Ro  \  J__ 
0  ?  'P  ' '  ?'a 

rt  G  L  rti 
o 
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This  is  the  normal  practice.  The  variations  in  transmitted  power,  P  and 
2 

Gq  are  found  using  a  delay  line  calibration  scheme  and  standard  point-target 
measurements,  respectively.  However,  these  approximations  of  constant-gain. 
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finite-beam,  constant-range,  constant-o0  inside  the  integral  may  not  be 
true  for  a  wide-beam,  an  antenna  pattern  with  high  sidelobes,  near-nadir  or 
near  grazing  angles  of  incidence. 

3.0  EFFECT  OF  ANTENNA  GAIN  AND  RANGE  VARIATION  INSIDE  THE  INTEGRAL 

As  a  first  attempt,  the  effects  of  a  more  realistic  antenna  gain  and 

4  o 

R  inside  the  integral  are  evaluated.  The  effect  of  a  variation  across  the 

illuminated  area  will  be  treated  in  Section  4.0. 

If  a  full  set  of  antenna  patterns  in  all  directions  is  available, 

these  numerical  values  can  be  entered  as  a  gain  function.  But  this  is  not 

usually  the  situation  and  some  kind  of  functional  approximation  of  the  antenna 

gain  is  inevitable. 

3. 1  Gaussian  Antenna  Gain 

If  the  antenna  gain  can  be  represented  as  a  Gaussian,  then 
G2  -  G2(ea,*a)  -  Go2  ^p(-2.773(ea-9p)2/Be23dB)exP(-2.773^a2/6a23dB)  (10) 

where : 

9p  =  antenna  boresight  pointing  angle  off  vertical  (in  x-z  plane) 

6  =  angle  off  vertical  in  elevation  plane  for  observation  direction 

d 

<j>  =  angle  off  x-z  plane  for  observation  direction 

d 

Be  =  two-way  3  dB  beamwidth  i >.  elevation  plane 

B  ...  =  two-way  3  dB  beamwidth  in  azimuth  plane 
a ,  iQo 

For  point  P(x,y)  on  the  ground  (see  Figure  3), 


0  *  tan  * (x/h) 


(11) 


(12) 


$  =  tan  (y//x2+h2) 
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This  coordinate  system  is  not  a  spherical  coordinate  system,  but  this  is 

used  because  the  antenna  patterns  are  available  in  these  azimuth  and  elevation 

planes. 


3.2  (sin  x/x;  Antenna  Gain 


One-way  antenna  gain  can  also  be  represented  by  (sin  x/x)  •  In  terms 


of  one-way  3  dB  beamwidth,  6 


3dB* 

s  in  (2.7831 6/BXjq)  2 


G(e)  ”  V  2.78310/6 


(13) 


3dB 


Therefore,  in  the  coordinate  system  shown  in  Figure  3»  two-way  antenna  gain 
can  be  represented  as 


,  ?  sinxfi  *  s  ' nX  .  * 

G  '  ‘o2<V’  V 


(14) 


where: 

X9  =  2.7831  (V8p^a,3dB 
X,  '  2-783'  *a/8a.3dB 

=  one-way  3  dB  beamwidth  in  elevation  plane 

e ,  juD 

3^  jjg  =  one-way  3  dB  beamwidth  in  azimuth  plane 
0  ,<!>  =  given  in  Eqs.  (11)  and  (12) 

3  3 


This  two-way  gain  function  given  by  Eq.  (14)  has  many  sidelobes  and  nulls, 
with  the  highest  sidelobes  26.5  dB  down  from  mainlobe  (see  Figure  4).  The 
location  of  the  first  sidelobes  will  be  about  1.6  away  from  the  center 


of  the  mainlobe. 


The  range  R  from  the  antenna  to  the  infinitesimal  area  element  dA  on 


the  ground  is 


R  =  /x2+y2+h2  (15) 

In  the  actual  HELOSCAT  III  system  (an  FM  radar)  the  range  to  the  target  is 
limited  by  the  bandpass  filter  in  the  final  stage.  Each  area  element  on  the 
ground  has  different  range  to  the  antenna  and  thereby  gives  a  different  IF 
frequency.  Therefore  a  bandpass  filter  with  finite  bandwidth  would  limit 
the  range  from  which  the  radar  receives  the  scattered  field.  In  the  narrow- 
beam  approximation,  the  bandpass  filter  plays  a  role  only  when  the  effective 
beamwidth  is  very  wide  and  the  incidence  angle  is  large  enough.  In  such  a 
case,  the  filter  cuts  the  area  of  the  ellipse  from  which  o°  is  calculated. 

The  transfer  function  of  the  bandpass  filter  shown  in  Figure  5  can  be 
represented  as  a  function  of  range  as  follows. 


f(R)  =  1, 

=  I  u  o 

-30(~  -(l+BV2f  )) 
=  10  Ro  0 


when  R  (1  -  B'/2f  )  <  R  <  R  (1  +  B'/2f  ) 
o  o  —  —  o  o 

when  R  <  R  (1  -  BV2f  ) 

—  o  o 

when  R  >  R  (1  +  B'/2f  )  (16) 

—  o  o 


where: 

R  is  the  distance  to  the  center  of  the  illuminated  area 
o 

B-*  is  the  bandwidth  of  the  filter  (Figure  5) » 
f  is  the  center  frequency  of  the  filter. 
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Theoretical ly  in  a  FM-CW  radar,  the  range  R  to  the  target  and  the  IF 
frequency  output  (F.^)  of  the  mixer  have  the  following  relationship. 


1»RAF  F 


Fif  " 


where: 

AF  =  RF  frequency  sweep  width 

F^  =  modulating  frequency 
c  =  velocity  of  light 

A  sample  of  the  transmitted  waveform  long  before  it  goes  out  of  the 
antenna  is  usually  used  as  the  local  oscillator  signal  to  be  mixed  with  the 
received  signal.  Hence  the  actual  IF  frequency  is  different  from  what  one 
expects  from  Eq.  (17),  especially  when  long  antenna  cables  are  used.  The 
discrepancy  in  range  should  be  obtained  experimentally  as  a  system  character¬ 
istic.  When  this  additional  range,  R  ,  is  not  considered  to  be  negligible, 

c 

Eq.  (17)  and  the  bandpass  filter  function  given  by  Eq.  (16)  should  be  modified 
accordingly  (see  Figure  5). 

1*(R+R  )AF  F 

r  cm  ,,0v 


f(R)  =  1,  when  R  (1  -B"/2f  )-(B'/2f  )R  <  R  <  (1+BV2*  )R  +(B'/2,  )R 
o  o  oc  —  —  too  toe 


,U  R+R 

(,-B'/2f0»>  «*«  *  <  0-»'/2fo)»o  -  (B-/2fo)Rc 

lU  o  c 


r*-rr\ 

=  I0'30(FTT'  "  0+B'/2fo)),  when  R  >  (l+BV2fQ)Ro  +  (BV2fQ)Rc  (19) 


Examining  Eq.  (19),  one  can  see  that  having  some  additional  range,  Rc»  has 
effectively  widened  the  bandwidth  of  the  filter. 


4-c 


$c  6* 

— 1 - 1 - 

Ro-t-Rc.  *+*<- 


7>(KK')  -r  Fr^uercu 


<J 


Figure  5s  Bandpass  Filter  Function 
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3.4  Integration  Limits 

In  this  section,  a°  will  be  assumed  to  be  constant  across  the  integration 
limits  and  evaluated  using  the  following  equation. 


a°  =  (4tt)3  f^r_j  _ 1 _ 

x2  pt  ||  G2(ea,<j>a)f (R)  ^ 


(20) 


In  view  of  the  narrow  beam  approximation  given  by  Eq.  (9),  one  can  rewrite  this 
equation  as 


II 


1 _ 

g2(ea,4>a)f(R) 


dA 


(21) 


where  G2(6  ,<fr  )  -  G  2g2(e  ,4  ). 

a  a  O  a  a 


The  integral  can  be  treated  as  a  weighted  area  [4]  compared  to  the  illuminated 
area  described  in  Section  2.0. 


w 


g  (ea,Of(R) 

_ q  a _ 


dA 


(22) 


The  integration  is  done  using  a  rectangular  coordinate  system.  There¬ 
fore  the  infinitesimal  area  dA  -  dx  dy  and  the  integration  limits  should  be  the 
whole  ground  plane  and  the  filter  function  will  effectively  set  the  limit. 

To  simplify  the  numerical  integration,  the  integration  limits  were  set  to 
be  the  points  where  the  antenna  gain  function  reduces  to  small  enough  numbers 
(about  -30  dB)  for  the  Gaussian  gain,  and  for  the  (sin  x/x)  gain  function 
the  limits  were  up  to  the  first  sidelobes.  In  this  case  the  integral  is 
also  evaluated  for  the  mainlobe  only  to  see  the  effect  of  sidelobes.  Figure 
6  shows  the  integration  limits  (square  area  marked  -30  dB  on  four  sides) 


together  with  the  elliptical  (shaded)  illuminated  area  normally  used  for  the 
narrow-beam  approximation.  Also,  the  equal-range  circles  which  effectively 
limit  the  area  by  the  filter  function  can  be  seen. 


The  weighted  area  now  becomes 
y+  x+ 


y-  x 


9  f0_(x,y),<j>  (x,y)]f[R(x,y)] 

a  a 

[R(x,y)/RQ]4 


dx  dy 


(23) 


where  g  ,  R  and  f  are  given  by  Eqs.  (10-12),  (14),  (15),  (19)  and 


x+  *  h  tan(6  +  1.6  6  ,JO)  for  Gaussian  gain 
p  —  e,3dB 

“  h  tan(6  +  3.2  8  ,  ._)  for  (sin  x/x)4,  up  to  first  sidelobes 

P  C  f jQ D 


and 


and 


h  tan(0p  +1.6  8e  j^g)  for  (sin  x/x)  ,  mainlobe  only 


y+  "  R^  tan(+  1.6  8  ,  .„)  for  Gaussian  gain 

—  o  —  a,3do 

■  R  tan(+  3.2  8  -JO)  ^or  (sin  x/x)4,  up  to  first  sidelobes 
o  —  a,  Jab 

*  R„  tan(+  1.6  8  ,.n)  for  (sin  x/x)4,  mainlobe  only 

o  —  a, 3dB  ' 


R  *  h/cos0 
o  p 


3.5  Results 

In  this  section,  the  possible  error  when  using  constant  range  and  con¬ 
stant  gain  across  the  effective  beamwidth  is  defined  as  follows. 

ERR0R1  =  A. /A  (24) 

1  w 

where: 
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Aw  Is  the  weighted  area  given  by  Eq.  (23)  and 

A.  is  the  area  of  the  ellipse  given  by  Eq.  (8) 

then 

ERROR]  (dB)  =  10  log (A. /A  )  (25) 

I  w 

and  this  ERROR!  (dB)  should  be  added  to  the  o°(dB)  calculated  using  a  narrow- 
beam  approximation. 

The  numerical  integration  was  done  using  Gauss'  formula  [5]  and  a  total 
of  1600  (40x40)  points  inside  the  integration  limit  shown  in  Figure  6.  Be¬ 
cause  the  number  of  points  evaluated  was  fixed  for  all  the  incidence  angles 
between  0°  to  70°  while  the  actual  integration  limits  in  x-  and  y-coordinate 
grew  tremendously  from  0°  to  70°,  the  accuracy  of  the  numerical  integration 
tends  to  get  worse  at  large  incidence  angles. 

Table  1  is  a  summary  of  ERR0R1  for  various  antenna  beamwidths  and  antenna 
gain  functions.  As  can  be  seen  in  the  table,  the  constant-range  and  constant- 
gain  across  the  effective  beamwidths  are  not  bad  approximations  to  use  up  to 
the  pointing  angle  of  60°,  except  when  the  two-way  beamwidths  are  as  large 
as  9°,  if  the  variations  of  a0  inside  the  integral  can  be  neglected  as  assumed 
here. 


4.0  EFFECT  OF  g°  VARIATION  INSIDE  THE  INTEGRAL 

Up  to  now,  c°  was  assumed  to  be  constant  within  the  limits  of  integration, 
so  ct°  was  outside  the  integral.  Actually,  o°  itself  is  a  function  of  the 
local  incidence  angle  and  therefore  a  function  of  x-  and  y-coordinates  of  the 
point  evaluated  (see  Figure  7). 
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TABLE  1 
ERR0R1 


Pointing 

Angle 

9p(°) 

Two-Way  3dB 

Beamwidths 

6e  =  9.0d" 
ea  =  9.o° 

=  6.3“ 

6  -  5.7- 

o 

"B  =  5.0 u 

B  =  5.0° 

3 

3e  =  274°” 

B  =  2.0° 

a 

Gaussian 

Gain 

0 

10 

20 

30 

40 

50 

60 

70 

0.055(dB) 

0.057 

0.060 

0.069 

0.095 

0.208 

0.591* 

0.366* 

0.025 (dB) 
0.025 

0.027 

0.031 

0.038 

0.063 

0.222 

0.616 

0.017(dB) 

0.018 

0.019 

0.021 

0.026 

0.035 

0.094 

0.588 

0.004(dB) 

0.004 

0.004 

0.005 

0.006 

0.008 

0.012 

0.033 

0 

0.059 

0.029 

0.022 

0.008 

4 

10 

0.060 

0.029 

0.022 

0.008 

i n  x j 

20 

0.067 

0.031 

0.023 

0.008 

30 

0.076 

0.037 

0.026 

0.009 

up  to 

40 

0.094 

0.046 

0.033 

0.010 

f  i  rst 

50 

0.228 

0.070 

0.044 

0.012 

sidelobes 

60 

0.583* 

0.206 

0.091 

0.021 

70 

0.484* 

0.484 

0.528 

0 

0.079 

0.049 

0.043 

1* 

10 

0.080 

0.050 

0.043 

0.030 

(^-^) 

20 

0.084 

0.052 

0.044 

0.030 

X 

30 

0.092 

0.056 

0.047 

0.030 

main  lobe 

40 

0.111 

0.063 

0.051 

0.031 

only 

50 

0.208 

0.080 

0.060 

0.033 

60 

0.561* 

0.213 

0.102 

0.038 

70 

0.340* 

0.582 

0.557 

0.052 

These  values  may  not  be  very  accurate  due  to  the  large  spread  between 
points  actually  evaluated. 


i 


.a-  ^ 
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o°  -  a°(0)  ■  a°(x,y) 

where  8  is  the  local  incidence  angle  and 


One  way  of  evaluating  o°  suggested  in  [3]  is  to  apply  Eq.  (21)  to  a 
set  of  measurements  to  get  a  rough  estimate  of  a°  variation  with  incidence 
angle,  then  use  this  function  in  the  integral  (Eq.  (1)),  and  continuously 
refine  the  o°  variation  function  to  get  better  estimate  of  a0. 

Variation  of  a°  with  incidence  angle  will  take  many  forms  depending  on 
the  kind  of  target.  In  many  cases,  including  sea  ice,  it  fits  quite  well 
with  the  exponential  function  used  here.  From  Eqs.  (1)  and  (21), 


and  o°(6)  for  some  particular  angle  of  incidence  can  be  evaluated  as. 


■5V 


g  f  [a,  (e)/o1 (e  )] 


(R/Ro) 
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M 


-b, (e-e  )  -b_  (e-e  ) 

Oj(0)/a°(0p)  =  e  p  and  a°(0)/a°(0p)  =  e  P 

and  ct°  variation  function  around  the  pointing  angle,  0p,  can  be  defined  as 

a  °!?(e)  “b  (0-0) 


A  W  / 

S  (0,0  )  =  ~ - 

n  P  a°(0  ) 


n  p' 


and  the  n-th  refined  o°  becomes 

4 

/i  \3  P  R 

o>„)  ■  •*¥-  ^X-T-)  - 5 - '■ - 

x  '  Go  ff  9  f  s„-l<9-9p> 


<R/Ro) 


Hence  the  integral  can  be  treated  as  another  weighted  area  similar  to  that 


of  Eq.  (23). 


g2  f  s  - 

A  4  1 - 0-idA 

ws  (r/rT 

o 


The  possible  error  when  the  a  variation  function  is  neglected  is  defined 


ERR0R2 (d B)  «  10  log(Aw/Awg)  (33) 

and  numerically  evaluated  (Table  2)  for  several  antenna  gain  functions,  beam- 

widths  and  a°  variation  across  the  integration  limits  (assuming  that  the 

variation  is  known  in  advance).  In  Table  2,  the  functions  of  e  and 

■0/9  k° 

e  *  were  derived  from  linear  regression  (in  dB  scale)  of  the  data  points 

obtained  from  sea  ice  using  the  narrow  beam  approximation  and  the  steep 

”0/5° 

variation  function  of  e  is  appropriate  to  some  open  water. 

In  the  next  section,  the  iteration  will  be  tried.  The  error  caused  by 
not  considering  steep  o°  variation  can  be  as  large  as  3.6  dB  for  normal 
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TABLE  2 
ERR0R2 


q°-variation  and  Two-way  3  dB  Beamwidths 


Pointing 

- - 

a  *  e 

~ "0  -0/12.9“ 

a  =  e 

7  -  3-e/5" 

0  _  -e/9.4u 

g  =  e 

Angle 

0  (°) 

Be  =  9.0° 

&  =  6.3° 
e 

Be  =  5.0° 

8e  =  2.4° 

p 

B.  =  9.0° 

8=  =  5.7° 

8,  =  5.0° 

8,  =  2.0° 

a 

a 

a 

a 

0 

3.647(dB) 

1 ,040(dB) 

2.140 (dB) 

0.532 (dB) 

10 

-0.484 

0.004 

-0.193 

-0.010 

Gaussian 

Gain 

20 

-0.994 

-0.060  . 

-0.314 

-0.020 

30 

-1.170 

-0.087 

-0.365 

-0.025 

40 

-1.301 

-0.108 

-0.401 

-0.029 

50 

-1.456 

-0.133 

-0.438 

-0.033 

60 

-1.386 

-0.190 

-0.505 

-0.038 

70 

-0.398 

-0.184 

-0.548 

-0.050 

0 

3.613 

1.026 

2.114 

0.526 

ij 

10 

-0.450 

0.008 

-0.184 

-0.010 

/Sin  X\ 

'  x  ' 

20 

-1.015 

-0.061 

-0.326 

-0.020 

30 

-1.330 

-0.089 

-0.377 

-0.025 

up  to 

40 

-1.405 

-0.111 

-0.417 

-0.029 

first 

50 

-1.359 

-0.122 

-0.437 

-0.033 

sidelobes 

60 

-1 .366 

-0.176 

-0.458 

-0.040 

70 

-0.633 

-0.205 

-0.555 

-0.044 

3.574 

1.012 

2.090 

0.517 

-0.466 

0.003 

-0.177 

-0.009 

-0.904 

-0.055 

-0.289 

-0.019 

A 

ma i n lobe 
only 

30 

-1.063 

-0.080 

-0.336 

-0.023 

40 

-1.177 

-0.099 

-0.369 

-0.027 

50 

-1.339 

-0.120 

-0.402 

-0.030 

60 

-1.388 

-0.177 

-0.459 

-0.035 

70 

-1.010 

-0.186 

-0.535 

-0.044 
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incidence  and  in  the  order  of  1  dB  for  other  incidence  angles  when  the  beam- 

widths  are  9°.  When  the  beamwidths  are  5°,  the  error  is  2.i  dB  for  normal 

incidence  and  less  than  0.5  dB  for  other  incidence  angles.  When  the  a0 

“0/12  9°  “0/9  4° 

variation  is  not  so  rapid  (e  or  e  ),  the  error  reduces  to 

1  dB  or  less  for  normal  incidence  and  becomes  very  small  for  other  incidence 
angles.  Recall  from  Section  3.5  that  the  numerical  integration  may  have 
some  error  for  large  incidence  angles.  Table  2  also  shows  that  the  (sin  x/x) 
gain  function  has  such  a  low  sidelobe  level  (-26.5  dB)  that  integrating  up 
to  first  sidelobes  did  not  differ  significantly  from  integrating  only  the 
mainlobe.  Furthermore,  the  (sin  x/x)  antenna  gain  function  and  the  Gaussian 
gain  function  did  not  show  very  much  difference. 

k. 1  I teration 

As  explained  in  the  previous  section,  iteration  is  tri*te*  with  r'-.?  1  data 
and  some  arbitrary  data  which  fits  well  with  the  exponential  function.  Be¬ 
cause  the  real  data  points  are  from  10°  to  70°,  only  these  points  are  evalu¬ 
ated  to  get  an  exponential  fit  and  then  iterated  several  times.  Table  3  is 
a  summary  of  the  result.  The  first  row  shows  the  o°  calculated  using  the 
narrow  beam  approximation  (0°^)  and  the  second  row  is  the  modified  a°  computed 
using  gain  function  and  range  variations  (o°^  +  ERR0R1 ) .  These  points  are 
regressed  to  get  an  exponential  fit,  and  this  exponential  function  is  included 
in  the  integration  to  give  the  third  row  (0°^  +  ERR0R1  +  ERR0R2).  From  these 
points  we  get  new  slope  and  therefore  new  ERR0R2,  and  for  the  fourth  row, 
o°  *  a0.  +  ERRORl  +  ERR0R2,  and  so  on. 


V 


-23- 
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TABLE  3 

DATA,  ORIGINAL  AND  UP  TO  4th  REFINED 


1.  13.6  GHz  Sea  Ice  Data  (ct°  =  e  ) 

»e  ■  2.«‘  sa  -  2.0“ 


ANGLE 

1C* 

'  20* 

30* 

40*  50’  60*  70* 

SLOPE  CORRELATION 

Go-U-<&  i 

2-70 

2.71 

2.7C 

2.70 

2.70 

2.60 

2.61 

2.59 

2.59 

2.59 

-1 .70 
-1.69 
-1.72 
-1 .72 
-1 .72 

-6.80  -10.00  -14.90  -26.CO 
-6.79  -9.99  -14.83  -25-96 
-6.82  -10.02  -14.92  -26.00 
-6.82  -10.02  -14.92  -26.00 
-6.82  -10.02  -U.92  -26.00 

-0.4621 

-0.4617 

-0.4622 

-0.4622 

-0.4622 

-0.9643 

-0.9644 

-0.9645 

-0.9645 

-0.9645 

ANGLE 

10" 

20’ 

30’ 

40*  50‘  60’  70' 

SLOPE  CORRELATION 

(m* 

“5i  «lc.lobe.$ 

2.70 

2.70 

2.69 

2.69 

2.69 

2.60  -1.70 

2.60  -1.70 

2.58  -1.72 

2.58  -1.72 

2.58  -1.72 

-6.80  -10.00  -14-90  -26.00 
-6.79  -9-99  -14.89  -25.97 
-6.82  -10.03  -14.93-26.02 
-6.82  -10.03  -14.93  -26.02 
-6.82  -10.03  -14.93  -26.02 

-0.4621 

-0.4618 

-0.4623 

-0.4623 

-0.4623 

-0.9643 

-0.9644 

-0.9645 

-0.9645 

-0.9645 

2.  4.8  GHz  Sea  Ice  Data  (o°  =  e_0/'2. 9“j 

Be  -  6.3“  B,  -  5.7" 


6-cu^fo^ 

(S-cmVv 

ANGLE  10*  20"  30*  40"  50’  60*  70° 

SLOPE  CORRELATION 

-2.70  -4-90  -10.50  -12.30  -14.80  -18.00  -24.00 
-2.67  -4.87  -10.47  -12.26  -14.74  -17.78  -23-33 
-2.67  -4.93  -10.55  -12.36  -14.87  -17.96  -23-56 
-2.67  -4.93  -10.55  -12.37  -14.87  -17-97  -23-56 
-2.67  -4-93  -10.55  -12.37  -14.87  -17-97  -23-56 

-C-3371  -C.9886 
-C.3293  -C.9893 
-C.3324  -0.99C1 
-0.3324  -C.99C1 
-0.3324  -0.9901 

cr ; , 

slacks 

ANGLE  10*  20°  30*  40°  50’  60°  70° 

SLOPE  CORRELATION 

-2.70  -4-90  -10.50  -12.30  -14-00  -18.00  -24. CO 
-2.67  -4.87  -10.46  -12.25  -14-73  -17-79  -23-52 
-2.66  -4.93  -10.55  -12.36  -14-35  -17-97  -23-72 
-2.66  -4.93  -10-55  -12-36  -14-85  -17-97  -23-72 
-2.66  -4-93  -10.55  -12.36  -14-85  -17.97  -23-72 

-0.3371  -0.9386 
-C.3309  -0.9894 
-0.3341  -0.9396 
-0.3341  -0.9396 
-C.3341  -0.9396 

3.  Arbitrary  Data  (o°  *  e’9^  fi  ■  6  =  9“) 


65oUU'»S  i 

ANGLE  10*  20*  30°  40“  50*  60*  70* 

SLOPE  CORRELATION 

-8.70  -17-00  -26.60  -33-OC  -42.00  -53.00  -60. CC 
-3.64  -16.94  -26.53  -32-91  -41-79  -52.41  -59-63 
-9- IT  -17. <30  -27.67  -34-17  -43-22  -53-73  -60.62 
-9-13  -17-93  -27.69  -34-20  -43-25  -53-80  -60.63 
-9-13  -17.93  -27-69  -54.20  -43-25  -53.80  -60.63 

-0.8618  -0.9987 
-0.8542  -0.9989 
-0.8636  -0.9989 
-0.8636  -0.9939 
-0.8636  -C.9989 

It,  T  VP 

up  te  fix i. 

'  “  1  ,1  J 

‘v'iA.it  r  * 

ANGLE  1C’  20*  30*  40“  50“  60*  70* 

SLOPS  CORRELATION 

-8.70  -17.00  -26.60  -33-OC  -42.00  -53-00  -60. CO 
-8.64  -16.93  -26.52  -32.91  -41.77  -52.42  -59-52 
-9-07  -17-92  -27-31  -54.27  -43-10  -53-76  -60.14 
-9.08  -17-93  -27.83  -34.29  -43.11  -53-77  -60.14 
-9- 08  -17-93  -27.93  -34-29  -43-11  -53.77  -60.14 
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For  all  the  cases  in  Table  3,  the  data  seems  to  converge  after  three 
iterations,  so  it  did  not  require  further  iteration.  Also  from  Table  1  and 
Table  2,  it  can  be  seen  that  the  ERRORl  and  ERR0R2  tend  to  cancel  out  each 
other  to  make  o°^  a  tolerable  estimation  at  least  for  the  incidence  angles 
of  10°  to  70°  and  when  the  beamwidths  are  not  very  wide. 


5.0  SUMMARY 

Three  kinds  of  so-called  illuminated  areas  are  treated  in  this  report: 

(1)  the  area  of  an  ellipse  with  constant  gain  and  constant 
range  to  the  antenna  for  all  the  cells  within  the 

el  1  ipse,  A. ,  (Eq.  (8)) . 

(2)  the  area  weighted  with  ^.itenna  gain  and  range,  Aw  (Eq.  (23)). 

(3)  the  area  weighted  with  antenna  gain  and  range  and  also 

a°  variation,  A  ,  (Eq.  (32)). 
ws 

The  possible  error  when  using  A.  instead  of  A^  was  defined  as  ERRORl 
and  summarized  in  Table  1.  When  the  o°  variation  with  incidence  angle  is 

small,  the  result  shows  that  the  constant  range  and  constant  gain  approxi¬ 

mation  gives  a  result  within  0.6  dB  or  less. 

Next,  a°  variation  is  included  and  the  possible  error  of  using  A 

w 

instead  of  A  was  defined  as  ERR0R2  and  listed  in  Table  2.  For  normal  inci- 
ws 

dence,  one  can  have  up  to  3.6  dB  error  when  the  beamwidth  is  9°  for  rapidly 

“ 0 /C°  •  Q 

varying  (e  3  )o  .  The  big  error  at  normal  incidence  is  always  present  and 
one  must  consider  this  problem  when  a0  for  normal  incidence  is  evaluated. 

The  error  will  be  bigger  than  the  values  obtained  here  if  the  real  a°  varia¬ 
tion  near  vertical  is  steeper  than  the  exponential  function,  which  is  often 
true  for  very  smooth  targets. 


The  more  realistic  00  values  corrected  with  these  estimated  errors 
were  further  refined  using  iterative  scheme.  Actually,  however,  the  iter¬ 
ation  did  not  show  any  further  improvement  because  it  converged  almost 
instantly  for  the  exponential  a0  variation  function  assumed  in  this  report. 
This  may  not  be  true  for  other  kinds  of  a0  variation  functions  and  can  be 
further  investigated  using  the  technique  described  in  this  report. 
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APPENDIX 
PROGRAM  LISTING 


IOC 

20C 

30C 

m 

60C 
70  C 

80  C . 

90C 
100C 
HOC 
J120C 
1  30C 
140C 
150C 
J60C 
170C 
180C 
1  90C 
200C 

21  gc 

220C 
230C 
.2  40C 
25CC 
260C 
270C 
280C 
290C 
300C 
31 OC 
320C 
330C 
340C 
350C 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
4908 
5008 
5108 
5208 
530 


5508 

5608 

5708 

58CC 

590C 


THIS  PROGRAM  IS  TO  EVALUATE  THE  POSSIBLE  ERROR  WHEN 
USING  THE  NARROW  BEAM  APPROXIMATION  IN  CALCULATING 
THE  SCATTERING  COEFFICIENT. 

ANTENNA  GAIN  F UN C T I  ON S / G AUS S I  AN  AND  (SIN(X)/X)**4, 

AND  REAL  RANGE  AND  TRAPEZOIDAL  8AND  PASS  FILTER  FUNCTION 
AND  THE  EXPONENTIAL  SIGMA  VARIATION  FUNCTION  IS  EVALUATED 
IN  THE  INTEGRAL  USING  GAUSS*  FORMULA.  40  POINT  BY  40  POINT 
TOTAL  1600  POINTS  ARE  EVALUATED  IN  THE  INTEGRATION  LIMIT. 

FUNCTIONS  REFERENCED 

GAIN1=GAUSSI AN  ANTENNA  GAIN 
GAIN2=(SIN<X)/X)T**4  GAIN 
BPF1=BANDPASS  FILTER  FUNCTION 
ASIGMA=EXPONENTIAL  SIGMA  VARIATION  FUNCTION 
R  EG  RES  =  L I  NEAR  REGRESSION  ROUTINE 
AVG=AVERAGI NG  ROUTINE 

ARE A= EFFECTIVE  AREA  CALCULATING  ROUTINE  FOR  THE 
NARROW  BEAM  APPROXIMATION  CASE 
VARIABLES  USED 

BE=TWO  WAT  3  DB  ELEVATION  BEAMWIDTHS 
B A« TWO  WAY  3  DB  AZIMUTH  BEAMWIDTHS 
A  NG  LE  =  PO I  NT  I NG  ANGLE 
T  HE0=ANGLE  IN  RADIANS 

W(40),X(40)=GAUSS"’  CONSTANTS  FOR  INTEGRATION 
R0=  RANGE  TO  THE  BEAM  CENTER 

UPX,UPY=UPPER  LIMITS  IN  X  AND  Y  COORDINATES 
LOWX,LOWY=LOWER  LIMITS 
H  =A  NT  £  NNA  HEIGHT 

ERROR1=ERROR  FOR  NOT  COUNTING  GAIN  AND  RANGE  VARIATION 

ERR0R2=ERR0R  FOR  NOT  COUNTING  SIGMA 

EAREA=EFFECTI VE  AREA  OF  THE  ELLIPSE 

WAREA=GAIN  AND  RANGE  WEIGHTED  AREA 

RAREA=SIGMA, GAIN, RANGE  WEIGHTED  ARES 

A L=  CONST  TO  SET  THE  LIMIT  OF  INTEGRATION 


COMMON  BE  ,B  A 
REAL  BE. BA 

REAL  ANGLE,THE0,PI,W<40)  ,X  (  4 0 ) , RO, UP  X. LO WX .UP Y, LO WY 

REAL  HALFX,HALFY,AVGX,AVGY,ISUM,JSUM,XX,YY,H,G,F 

REAL  I SU Ml, JSUM1.PROD1.ERROR1 (7). ERR  OR  2(7), RAREA 

R EALP ROD, RANGE, INTGRL, ERROR, ERR D3 

REAL  WAREA,EAREA,AREA,ABE,ABA,EB£,EBA 

REAL  FILTER, BPF0.BPF1 

REAL  B CONS T,C CONST, SI GMA , A L , A S I GMA 

REAL  SIGMA 3(7), SIGMA1 (7) 

INTEGER  I,II,J,K 
DATA  H/3.6/ 

DATA  <W(I) ,1  =  1,20)/. 3 7750594,. 077339 81,.G 7 61  1036,. 374 72 31  6, 

.3728865 8,. 07061 1 64 ,. 36 791 2 3 4 ,. 36 4 80401 , 
.06 1306 24,. 0574 39 76,. 0532 2784,. 34 869560. 
.04387090,. 03 8782 16,. 0334601 9, .32793700, 
.0222 45 8 4,. 01 64 2 105,. 01 04 9828,. 004 52 127/ 


DATA 


(  X(  I), 1  =  1, 20)/. X3  387  72  41, .  11  63  54  07, ,  1  9269758. .26  81  5218 
.34199409, .41 3779 20..4 6 307 560,. 5 4 946 71 2 
.61 25 5388,. 6 71 9 5 068,. 7 273 1825,. 77830565 


* 


.8246 122 3,. 86 59 59 50,. 90209880,. 93281 280, 
.95 79 1651,. 97725 994,. 9907 26 2 3,. 99823 770/ 


SYSTEM  CONST  FOR  INTERNAL  RANGE 
OATA  R CONST/2. 6/ 


61  OC 
620C 
630 

fs8c 

660C 

670 

..680 

690C 

700 

710 

720C 

730 

740 


CONSTANT  TO  CONVERT  SLOPE 
B8B=10.*AL0G10(EXP(1«.)) 
PI*4.*ATAN(1.0) 


IN  OB  SCALE  INTO  EXPONENTIAL  SLOPE 


TWO  WAY  3 
ABE*  9. 0 
ABA*  9. 0 


B  E*A  8  E  *P  I / 1 80. 
BA* ABA *P 1/180. 
TWO  WAY  EFFECTIVE 
E BE  =  1 . 201  * A  BE 
EBA* 1.201  * AB A 


OB  BEAMWIDTHS  IN  DEGREES 


6  EAMW 1 0  THS  IN  DEGREES 


699 


_  76. 

"  770 
780 
790 
800 
'810 
820 
830  698 

_ 84  08  70 

850 
860 
870C 
880C 
890C 
900 
910 
920 
930 


PARAMETERS 

-42./-53W-60./ 


697 


INPUT  SIGMAO  AND  OUTPUT  THESE  INPUT 
DATA  SIGMAO/ “8. 7,-1  7.., -2  6.  6  3 3  .  , 

N=7 

CALL  REGRES(SIGMAO,N, SLOPE, CORR) 

WRI T  E  (  6,  699) 

FORMAT  (/IX, 'SIGMAO  DATA,  ORIGINAL  AND  UP  TO  4 T H  REFINED') 
WRITEC6/698) 

F0RMAT(/1X, 'ANGLE  10  20 

SLOPE  CORRELATION') 

WRITE (6,697)  (S IS*AO< 1 > ,  1*1 ,7) , SLOPE ,CORR 
F0RMAT(4X,7F7.2,1X,F7.4,2X,F7.4) 


40 


50 


60 


94QC 
950C 
960C 
970 
980 
990C 
1000 
1010 
1020 
1030C 
1040  C 
1050 

I  060  C 
1 070  C 
1080 
1090 

1H8 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

I I  90  C 
1200 


ASSIGN  INTEGRATION 
DO  100  1=21,40 
W ( I)=W(I-20) 

X  (  I)=-X(I-20) 
100  CONTINUE 


CONSTS 


DO  UNT IL  4  TH  ITERATION 

DO  89  11=1,4 
SLOPP=—BBB/ SLOPE 

DO  90  K*1 , 7 
ANGLE=FL0AT(K)*10. 

THE0=PI*ANGLE/180.. 

RO*H/ COS (THEO) 

LIMITS  UP  TO  FIRST  SIDELOBES  OF  ( S I N ( X) / X ) *  *  4 
AL=PI*1. 414/2. 7831 
ALU*A  L 

IF(K  .EQ.  6  .OR.  <  .£Q.  7)  ALU=1.5 
UPX=H*(TAN(THEO+ALU*BE)) 

LOWX*H*<  TAN( THEO-AL*BE) ) 

UPY=RO*TAN( AL*BA) 

LOW Y*-UPY 

HAL  F  X  =  (UPX-LOWX> /2. 

HAL  FY  =  ( UPY-LOWY) /2. 

AVGX*  (UPX*L0WX)/2. 

AVG  Y*  (UPY  +  LOWY) /2. 


ISUM1 =0 


1210 

1220 

1230 

1IS8 


I$UM*0. 

DO  101  1*1/20 

YY*X< I ) *HALF Y+AVGY 
J  SUM*0, 

J  SUH1 =0. 


1 260  C 

1270  DO  102  J  *  1  /AO 

_ 1280 . .  .  XX=X(J ) *HALFX+AVGX 

1290  CALL  GAIN1 (XX/YY/H/ THEO/G) 

1300  RANGE=SQRT(XX**2+YY**2+H**2) 

1310  FILTER=BPF1<RANG£/R0/RC0NST) 

_ 1320  RANGE4*RANGE**4 

1330  SIG MA*A SIGMA (XX/YY/H/ANGLE/SLOPP) 

1340  PROD=G*F I L T E R * S I GMA 7 R ANG E 4 

1  350  C 

_ 1  360  PR0D1=G*FILTER/RANGE4 

1370  JSUM1*JSUM1*W(JJ *PR0D1 

1380  JSUM=JSUM*W(J)*PROD 

1390  1  02  CONTINUE 

_  1 400  C 

'1410  I  SUM=ISUM  +  W< I) *JSUM 

1420  ISUMl=ISUMl+W<I)*JSUM1 

1430  101  CONTINUE 

1  440  C 

'1450  iNTGRL=2„*HALFX*HALFY*ISUMl 

1  4  60  C 

1470C  WEIGHTED  AREA  USING  3 DB  BEAMWIDTHS  WITHOUT  SIGMA  VARIATION 
_1 480  WAREA  =  R0**4*INTGRL 

1 4  90  C 

1 500  C  WEIGHTED  AREA  INCLUDING  SIGMA  VARIATION 
1510  RAREA=2.*HALFX*HALFY*ISUM*(R0**4) 

1  520  C 

1 530  C  EFFECTIVE  AREA  USING  EFFECTIVE  BEAMWIDTHS 

1540C  AND  CONST  GAIN  AND  CONST  RANGE 

1550  EAREA=AREA(EBE/EBA/ANGLE/H/RCONST) 

1.560  C 

1570  ERR0R1(K)*1O.*ALOG10(WAREA/EAREA) 

1580  ENROR200  =  10.*AL0310(RAREA/WAREA> 

1  590  C 

1 600  C  WR I T  E<  6 / 999 )  ANGLE/ ERRORl /ERR0R2 

1610C  999  FOR  MAT (1X/*ANGLE** /f5.1/2X/,ERR0R1='/F6.3/'  DB*/2X/ 

1 620  C&  •ERR0R2*’ /F6.3/'  DB’) 

1 630  C 

1 640  C  WRITEC6/997)  EARE A/ WAR EA/ RARE  A 

1 6  5  0  C  997  FORMAT (1 2X/ 'EAREA** / F9.3/ 2X/ ' WARE A  =  '/F 9.3/ 2X/' RARE A=* /F9.35 
1660  C 

1670  90  CONTINUE 

1680  C 

1690  1 F  (  1 1  .NE.  1)  GC  TO  112 

1700  00  111  1*1/7 

1710  SIGMA1 (I)«SIGMAO<I)-ERROR1 (I) 

1720  111  CONTINUE 

1 730  C 

1740  CALL  REGRESCSIGMAI /N/SLCPE/CORR) 

1750  WRI TE  <6 / 697 )  < S I G M A  1 ( I ) , I s 1 , 7  ) / SL 0 P E / CO R R 

1 7o0  GO  TO  114 

1770  1  12  CONTINUE 

1  780  C 

1790  00  113  1*1/7 

1800  SIGMAO(1)*SIGM41CI)-ERROR2(I) 


113  CONTINUE 


3118  c 

1830 

1840 

1850 

1860 

1870 

-3118  c 

1900  C 
1910 
_1 920  C 
1  930  C 
1  940  C 
1  950  C 
I960 
1970 
1980 
1990 
2000  C 
2010 
2020 

1818 
2050 
2060 
2070 
2080 
2090  C 

i!?8c 

21  20  C 
21 30  C 
2140C 
21  50  C 


CALL  REGRES<SIGMAD*N*SLOPE*CORR> 

WRI  TE  (6*697)  < S I GM AO < I > / I  =  1  *  7 > *SLOPE / CO R R 
114  CONTINUE 
89  CONTINUE 
STOP 
END 


2170 
2180 
2190 
22C0C 
2210  C 
2220 

Illgc 


2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

23  70  C 

2380 

2390 

2400 


SUBROUTINE  GAIN1 (X*Y*HH*TO#GG) 

GAUSSIAN  ANTENNA  GAIN  FUNCTION 
WITH  TWO  WAY  BEAMWI0THS=8E/BA 

COMMON  BE/BA 

REAL  BE*  BA 

REAL  XsY/HHsTO/GG 

REAL  THEA,PHIA*G1*G2sAA 

THEA«ATAN(X/HH) 

PHIA«ATAN(r/<SQRT<X*X+HH*HH)>> 

AA«4.*AL0G<2.) 

G1«EXP(-AA*< (THEA-TO) **2)/ (BE** 2)  ) 
G2*EXP<-AA*<PHIA**2J7<3A**2>  > 

GG=  G1  *G  2 
RET  URN 
END 

SUBROUTINE  G A I N2 < X  * Y *HH *T 0*6G > 

(  S  IN ( X) /X ) *  *4  TWO  WAY  ANTENNA  GAIN. PATTERN 
WITH  TWO  WAY  BE AMW I D TH S »B E * B A 

COMMON  BE/BA 

REAL  BEsBA 

REAL  X*  Y  *  HH  *  TO  *  C  G 

REAL  THEA,PH1A,G1*G2*BE0NE*BA0NE/A1*A2,G11,G22 

ONE  WAY  3 DB  BEAMWIDTHS 
BE0NEeBE*SQRT (2.) 

BA0NE=6A*SGRT (2.) 

THEA=ATAN(X/HH) 

PHI  A= AT  AN <  Y/ ( S Q R T < X  * X  +  HH*H H ) ) ) 

A1»2.78  31  *<THEA-T0)/BE0NE 
A2*2. 7831 ‘PHIA/BAONE 
G1*ABS< CSIN< A1 ) } / A  1 } 

IFCGI.LE.  (1.0E-13))  GO  TO  400 
G1 1  *G 1  *  *  4 

G2*ABS< (SI N(A2) ) / A  2  > 

I F  (  G  2  .LE.  (1.0E-10))  30  TO  400 
G22  *G  2»  *  4 
GG*  G 1 1*  G22 
RETURN 


4G0  CONTINUE 
GG*  0. 

RET  URN 


2410 

Itl8c 

2440  C 
24  50  C 
2460C 
2470 
2480 


2510 
2520 
2530 
2540 
2550 
2560 
2570 
2580 
2590 
2600 
2610 
2620 
2630  C  C1 
2640 
2650  C 
2660  C 
2670  C 
2680 
2690 
2700 
27 1 0  C 
2720 
27  30 
2740 
2750  C 
2760 
2770 

f?!8  8< 
2800 
2810 
2820 
2830C  C ’ 
2840 
2850  C 
2860  C 
2870  C 
2880  C 


2920 

2930 

2940 

2950 

2960 

2970 

2980  C 

2990 

30C0C 


REAL  FUNCTION  B  PF 1 (R/R0/C0NST) 

TRAPEZOIDAL  BANDPASS  FILTER  FUNCTION 
REAL  R/RO, CONST 

REAL  FIL,BW1,CF,L0WR,HIGHR,RR0,HF,LF 
DATA  BW1,CF/12.C,50..0/ 

FIL*BW1 / <2. *CF> 

HF *  1 . +F  IL 
LF=1.-FIL 

LOWR=RO*LF-FIL*CONST 
HI6HR=R0*HF+ FI L* CONST 
RRO=(R+CONST)/(RO+CONSTO 
IF ( R.GT.HIGHR  .AND.  R R 0 . L E . < 1 . 2 *H F > )  BPF 1 =1 0. ** (“30. * ( RRO-HF ) ) 
IF (RRO.GT . (1 . 2*H  F  )  5  BP  F 1 =0  . 

IF< R. LE.HIGHR  .AND.  R.GE.LOWR)  8PF1=1. 

IF(R.LT.LOWR)  BPf 1 *1 0. *  * < 1 4 5 . * < RRO- LF) /4 . ) 

RETURN 
END 

REAL  FUNCTION  BPFO<R,RO/CONST) 

RECTANGULAR  BANDPASS  FILTER  FUNCTION 

REAL  R/RO/CONST, FIL 
REAL  BW/CF /LOWR/HIGHR 
DATA  BW/CF/13. 5/53.0/ 

FIL=8W/ <2.*CF> 

LOWRsRO*  <1 .-FIL)-CONST*FIL 
HIGHR=R0*(1 .+F IL)+CONST*FIL 

IF ( R. GE. LOWR  .AND.  R.LE.H1GHR>  GOTO  803 
BPF  0*0. 

RETURN 
)0  CONTINUE 
BPF  0=1  • 0 
RETURN 
END 

k*******  ****** 

REAL  FUNCTION  AS  I GMA < X, Y, HH, ANGL, B ) 

THIS  FUNCTION  IS  THE  SIGNA  VARIATION  FUNCTION  FOR  POINT 
P(X,Y)  WITH  POINTING  ANGlE=ANGL  AND  SLOPE=B. 


REAL  X, Y/HH, ANGL/B 
REAL  THETA/ATHETA/PI 
PI*4.  *AT AN (1 .) 

THE  T  A«AT  AN ( SQRT(X*X*Y*Y)/HH> 

ATHETA*THE  TA*1 80. /PI 

AS  I GMA*EXP(-( A  THE  T A -ANGL) / 6) 

RETURN 

END 


SUB  ROUT  I NE  REGRES<Y,NN, SLOP/ COR ) 


3010  C  THIS  SUBROUTINE  CALCULATES  THE  SLOPE  AND  THE  CORRELATION 
3020C  COEFFICIENT  OF  AN  ARRAY  Y(NN).  THE  X-VARIABLE  IS  THE 
3030  C  INCIDENCE  ANGLE  FRO*!  10  TO  10*NN  DEGREES. 

2040  C 

3050  INTEGER  NN 

3060  REAL  Y(NN), SLOP, SLOPE, COR, X(10),SY(10),SX(10) 

3070  REAL  SXY(100,AVG,AVX,AVY,9BB 

3080  C 

3090  DO  10  I  * 1 , NN 

3100  X< I )=FLOAT< I) *1 0. 

3110  10  CONTINUE 

3120  NNN*NN 

3130  AVX  =AVG  <  X,NNN) 

3140  AVY=AVG(  Y/NNN) 

3150  DO  11  1*1, NNN 

3160  SY(I)*(Y(I>-AVY)**2 

3170  SX(  I)*(XC  D-AVX)  **2 

3180  SXY (I)*(X(I)-AVX5*(Y(I)  —  AVY) 

3190  11  CONTINUE 

3200  SLOP=AVG  <SXY,NNNI/AVG<SX,NNN> 

3210  COR=SLOP*SQRT< AVG(SX,NNN) /AVG(SY,NNN)> 

3220  RETURN 

2230  END 


- 

5250 
3260  C 
32 70  C 
3280  C 
3290  C 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 

*381 c 

34  10  C 
3420  C 
3430  C 
3440  C 
34  50  C 
3460  C 

34  70  C 
3480  C 

3*81! 

3510  C 
3520  C 

35  30  C 
3540 
3550 
3560 
3570 
35  80  C 


************ 

REAL  FUNCTION  AVG ( A, N) 

THIS  FUNCTION  CALCULATES  THE  AVERAGE  OF  AN  ARRAY  WITH  N 
ELEMENTS. 

INTEGER  N 
REAL  A ( N ) , S  UM 
SUM  *0  • 

DO  12  I*1,N 
S  UM*S  UM  +  A ( I ) 

12  CONTINUE 
AVG=SUM/FLOAT(N) 

RET  URN 
END 

*********  **** 

REAL  FUNCTION  A  R  E  A  (  E  LB,  A  Z  3  ,  A  NG  L  ,H  ,  R  CONS  T  ) 

THIS  FUNCTION  AREA  CALCULATES  THE  EXACT  ILLUMINATED  AREA 
FOR  AN  FM-CW  RADAR.  CALCULATION  INCLUDES  AREA  LIMITING 
BY  THE  BAND  PASS  FILTER  . 

INPUT  PARAMETERS 

E  LB  :  E  L  EV  A  T I  ON  BEAMWIDTH  IN  DEGREES 
AZBiAZIMUTH  BEAMWIDTH  IN  DEGREES 
ANGL: INCIDENCE  ANGLE  IN  DEGREE 
HiHElGHT  OF  THE  ANTENNA 
RCONST {SYSTEM  CONST  FOR  INTERNAL  RANGE 

REAL  PI,H 

REAL  AN GL,ELB,AZB, RCONST, ELSR,A23R 
REAL  THE,FC,BW,FIL,RE,RA 

REAL  R,C,ALPA1,AL»A2,GAM,A,B,K,L1,L2,X1,X2 


3590  C  FUNCTIONS  REFERENCED 
3600  C 


u 


3620  C 
3630 
3640 
3650  C 
3660  C 
3670  C 
3680 
3690 
3700  C 
3710 
3720  C 
3730 
3740 
37  50 
.3.760 
3770 
3780 
3790 
3800 
3810 
3820 
3830 
.2 8.40  „ 
3850 
3860 
3870 
3880 
3890 
3900 
3910 
3920 
3930  C 
3940  C 
3950  C 
3960  C 
3970 
3980 
3990 
4000  C 
4010 
4020  C 
4030 
4040 
4050  C 
4060 
4070  C 
4080 
4090 
41  DOC 
4110 
4  1  20  C 
4130 
4140 
41  50  C 
4160 
41  70  C 
4180 
41  90 
4  200  C 


REAL  LAREA1 ,LARE A2 

DATA  FC/BW/50.0,13.5/ 

PI=3. 1415926 

CONVERT  BEAM WIDTHS  TO  RADIANS 

ELBR=PI*ELB/180. 

AZ3R=PI*AZB/180. 

I F (  IF  IX(ANGL)  . E  Q.  0)  GO  TO  301 

THE=ANGL*PI/180. 

R  =  H /C OS (THE) 

C  =  H*T AN<  THE) 

FIL=BW/ (2.*FC> 

ALPA1=R*(1.-FIL)-RC0NST*FIL 
ALPA2=R*<1„+FIL)+RC0NST*FIL 
GAM=H*(TAN(THE-ELBR/2.) ) 

RA»  2  .  *R*(TAN(AZBR/2.,)) 

RE=H*(TAN(THE+ELBR/2.))-GAM 
A  =  R  E/2. 

B=R  A/2. 

K*G AM+A 
L1*C-GAM 
L2=  RE— L 1 

IF<  ALPA1-H)  120/120,122 
120  XI  *  C 

GO  TO  124 

122  X 1=C-SQRT (ALPA1 *ALPA1 -H*H> 

124  CONTINUE 

X2*SQRT ( ALPA2**2-H*H)-C 

CHECK  IF  THERE  IS  ANY  BEAM  LIMITING  DUE  TO  NARROW  BPF 
AND  IF  ANY/  GO  TO  APPROPRIATE  CASES 

1 F ( XI -Li)  200/201/201 

200  IFCX2-L2)  205/203/203 

201  IFCX2-L2)  204,202/202 

202  CONTINUE 

CASE  1/  NO  LIMITING  OF  THE  BEAM 
AREA  =PI  *A*8 
RETURN 

203  CONTINUE 

CASE  2/  NEAR  SIDE  OF  THE  BEAM  IS  LIMITED  BY  THE  FILTER 
AREA=PI*A*B-LAREA1 (A/B/C/K/Xl  ) 

RETURN 

204  CONTINUE 

CASE  3/  FAR  SIDE  OF  THE  BEAM  IS  LIMITED  BY  THE  FILTER 
AREA=LAREA2<A/B/C/K/X2) 

RETURN 

205  CONTINUE 

CASE  4/  BOTH  SIDES  OF  THE  BEAM  IS  LIMITED  BY  THE  FILTER 
AREA=LAREA2<A/B/C/K/X2>-LAREA1(A/B/C/K/X1) 

RET  URN 


4210 
4220 
4230 
4240 
42  50  C 
4260  C 
4270 
.4280  C 
4 290  C 
4300  C 
4310  C 
4320 
4330  C 
4340 
4350 
4360 
4370 
4380 
4390 
4400 
4410 
4420 
44  30 
4440  C 
44  50  C 
4460 
44  70  C 
4  4  80  C 

44  90  C 
4500 
4510 
4520 
4530 
4540  C 

4 5  50  C 
4560 

4  5  70  C 
4580  C 
45  90  C 
4600 
4610 
4620 
4630 
4640 
4650 
4660  C 
4670  C 
4680 
4690  C 
4700  C 
4  7 1 0  C 
4720C 
4730 
4740 
4750 
4760 
4770 
47  80  C 
4790  C 
4800  C 


301  CONTINUE 

ARE  A  =  PI *H*H*TAN<ELBR/2. )  *  T  ANCAZBR/2.) 

RETURN 

END 


REAL  FUNCTION  I N TS EC < A1 , A2 / A  3 / A4) 

THIS  FUNCTION  CALCULATES  THE  INTERSECTION  Of  THE  ELLIPSE 
<X-K) **2/A**2*Y**2/3**2=1  AND  THE  CIRCLE  X  * *2 ♦ Y * *2 *A 4  *  *  2 

REAL  A1 * A2/A3sA4 
LOCAL  VARIABLES 

REAL  B1 sB2/B3sB4s35#B6 
B1* ( A2+A1 ) *< A2-A1) 

B2s(A4+A2)*(A4-A2) 

B3= <A 1*A2*A3)**2 
B4s a2*A2*A3 
B5=  A 1 *A 1 

B6=SQRT  <B3-B5*B1 *32) 

1NTSEC=(B4-B6)/31 

RETURN 

END 


********** 


******** 


REAL  FUNCTION  MISCCTHETA) 


THIS  FUNCTION  CALCULATES  < 2  *  THE T A-S I N < 2*T H E T A) > /4 


REAL  THETA 

MISC= (2. *THETA-S IN(2.*THETA> )/4. 

RET  URN 

END 


*********** 


******** 


REAL  FUNCTION  T  H  EX  <  MG*  MP) 


THIS  FUNCTION  CALCULATES  A R C T AN ( S3 R T ( MG *  * 2-M P *  *  2) / M P ) 


REAL  MG/MP 
REAL  MGP 

MGP*$QRT  <MG*MG-MP«MP> 

THEX=ATAN(MGP/MP) 

RETURN 

END 


************************************************ 
REAL  FUNCTION  L  A  R  E  Al  <  A  A  *  B  3  *  C  C /•  K  K*  X  X  1  ) 


THIS  FUNCTION  CALCULATES  THE  NEAR-SIDE-BEAM  LIMITED  AREA 
WHICH  SHOULD  BE  SUBTRACTED  FROM  AREA  CALCULATION 

REAL  AA*B6*CC*Kfc*XXl 

REAL  TTH1*IINT1*K<MP*TTHE2/IINT2*GG*PPSI 
REAL  THE  X* Mi  sc* I  NTS  EC 
GG=  C  C  -X  X  1 

PPS  I*INTSEC<AA*6B*MC*GG> 

INTEGRATE  ALONG  THE  CIRCLE 


■4810  TTHl=THEX(GG*PPSI> 

4820  IINT1=GG*GG*MISC(TTH1) 

4830  C 

4840C  INTEGRATE  ALONG  THE  ELLIPSE 
4850  C 

4860  KKMP=KK-PPSI 

4870  TTHE2=THEX(AA/KKMP) 

.4880  .  1IN T2=AA*8B *MI S C < TTH E2) 

4890  C 

4900  LAREA1*2«*(IINT1+JINT2) 

4910  RETURN 

_4920  END  _ 

4  930C 

4940C  ********************************************** 

4950  REAL  FUNCTION  L A  RE A2 < A A  *  3 3  * C C  /  K K, X X 2  ) 

. 4960  C 

4970  C  THIS  FUNCTION  CALCULATES  THE  AREA  FOR  THE  CASE  WHEN 
4980  C  THE  FAR  SIDE  OF  THE  BEAM  IS  LIMITED  BY  THE  FILTER 
4  990  C 

5000  REAL  AA*BB*CC*KK*XX2 

5010  REAL  G*PSI *THEl *THE3/ INTI ,INT2/PI*KMP 

5020  REAL  THEX,MISC*INTSEC 

5030  DATA  PI/3.1415926/ 

5040C 

5050  G*C  C  +  XX2 

5060  PS  IaI NT  SEC (AA*BB*KK*G) 

5070  C 

5080  C  INTEGRAL  FROM  PSI  TO  G  THROUGH  THE  CIRCLE 
5090C 

5100  THE1=THEX(G/PSI) 

5110  INT1*G*G*MISC<THE1 ) 

51  20  C 

5130C  INTEGRAL  FROM  GAMMA  TO  PSI  THROUGHTHE  ELLIPSE 
5140C 

5150  IF(KK-PSI)  300*301 *302 

51  60  300  KMP ®P S I  —  K< 

5170  THE3=THEX( AA*KMP) 

5180  I NT  2*AA*B3*{.5*PI“MI SC  <  THE  3) ) 

5190  GO  TO  303 

5200  301  INT2=PI*AA*BB/4. 

5210  GO  TO  303 

5220  302  THE  3= TH E X ( A  A *K < -PS  I ) 

5230  INT2=AA*BB*MISC(THE3) 

5240  C 

5250  303  LAREA2=2.*<INT1+INT2) 

5260  RETURN 

5270  END 


END 

FILMED 


